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ON THE FREE COOLING OF A HOMOGENEOUS SPHERE, OF INITIAL 

UNIFORM TEMPERATURE, IN A MEDIUM WHICH MAINTAINS A 

CONSTANT SURFACE TEMPERATURE. 

By Mr. R. S. Woodward, Washington, D. C. 

1. The problem of the cooling of a homogeneous sphere initially heated to 
a uniform temperature divides itself naturally into two cases, namely : first, that 
in which the dissipation of heat at the surface of the .sphere goes on independ- 
ently of the surrounding medium, and second, that in which the dissipation is 
conditioned by the surrounding medium. For brevity we have applied the phrase 
" free cooling" to designate the first case. Both cases have been stated by Fourier 
and subsequent writers, but their investigations are, so far as we know, confined 
almost wholly to the second and more complex case. Moreover, they have, for 
the most part, contemplated the subject from the standpoint of the pure analyst; 
they have been too much occupied with the difficulties of their splendid analysis 
to give much heed to the needs of the computer. 

The most interesting example for application of the theory of a cooling 
sphere is presented by the earth. The hypothesis most commonly entertained 
is, that the whole mass of the earth was at a certain epoch heated to a uniform 
temperature, and is now slowly cooling by conduction without sensibly heating 
surrounding space. With the ultimate object in view of tracing out the conse- 
quences of this hypothesis, we have sought to express the solutions of the case 
of free cooling and that of conditioned cooling in such terms that the computer 
can, without undue labor, assign the temperature at any point within the sphere 
for any value of the time. 

In the following pages the case of the free cooling subject to the restriction 
of a constant coefficient of diffiision is alone considered. How closely the 
conditions of this restricted case accord with those actually presented by the 
earth is a question which requires examination in its proper place. The present 
enquiry, however, is directed solely to the mathematical treatment of the prob- 
lem, and the conditions are tacitly assumed to apply to the earth. 

2. The law of cooling of a homogeneous sphere, initially heated to a uni- 
form temperature and cooling in a medium which maintains a constant surface 
temperature, must evidently be such as to give the same temperature for all 
points equidistant from the centre of the sphere. 

Let Uf, be the initial uniform excess of the temperature of the sphere over 
that of the surrounding medium. Let u be the excess of the temperature of any 
shell of radius r and thickness dr over the temperature of the surrounding me- 
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dium at any time t, after the initial epoch. Denote the radius of the external 
surface of the sphere by r„ ; and let the coefficient of diffusion which is here 
assumed to be constant, be symbolized by c?. Then the function of ?/o, u, r„, r, t, 
and (f which expresses the law of cooling must be such as to satisfy the follow- 
ing partial differential equation : — 

iM — ^2f!M (,\ 

If we write the functional relation which is to be discovered in the form 

ru ^f{u^, Vfi, r, t, a^), 

we must also have r« = o for r = r^, (2) 

ru = for r = o, (3) 

ru = rUf) for / = o, (4) 

ru ^=0 for (= cc . (5) 

Without going through the steps essential to exclude other functions, we 
observe that conditions (i) and (2) are satisfied by the expression 

Qe ~ "' (»^ /'■oV' sin n:zr I r,, 

C„ being any constant, and n any integer. Hence, since every term of this form 
will satisfy (i) and (2), we may write 

«/ =" I" \e - ''^ ("'"/ '■0) '' sin «;r r / r„ . (6) 

« = I 

When (=: o, (6) becomes in accordance with (4), 

« = CO 

ru„ = 2 Cn sin JtTz rj r„. (7) 

«^ I 

To fulfill this condition we must have 

'•0 ''0 

u^ I ;• sin nn ~ dr = C„ I sin^ nrr - dr. (8) 

o " o " 

r 
If we put d ^ TT ~ , 

''0 

(8) becomes '^ jO sin nd dd = C„ J'sin^ nd dO, 

o o 

, ^ 2r,. «„ cos KTt 

whence Q = -^^— ^ • 



7: n 

2r, u, (-1) " + ' 
n ' n 



(9) 
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Substituting this value of C^ in (6), we find the following equation, which ex- 
presses the law of cooling sought, in its simplest form : — 



2" 
« = I 



ru =z " " 2 ^ '- e \ I "I sin ^Tt — 



+ e ^10' sin TT- 



(lo) 









" > 



+ ^^ '-^ ' "' Sin 3;r- 



The equation (lo) satisfies all of the conditions (i) to (s), because each term 
satisfies them. 

When ^ = o. (lo) gives for the temperature at the centre of the sphere, 

2r„ «o o 2^0 ?/„ d (rti) 

Tt ' o TT ' dr ' 

Evaluating this expression, we find for the temperature at the centre of the sphere. 



When ^ = o this gives u = it„, as may be readily seen by recurring to (to). 
When ^ > o it appears that u <,u„. 

3. The series (lo) and (11) express all the circumstances of the cooling 
sphere under the assumed conditions. They are rapidly converging series when 
«^(;r/;'„)^^ is not much less than unity, since the negative exponents in the succes- 
sive terms increase as the squares of the natural numbers. But in the most im- 
portant application of this theory, namely, to the earth, the value of ^^(jr/;-,,)^/ is 
very small unless the time, /, is very large. Thus, for the earth, if the British 
foot and the year are the space and time units, we have in round numbers 

a^ = 400,* 
r„ ^ 2 1 000 000 feet. 
Hence, in order to make a^iit jr^^t unity, we must have 

(2 I OOP ooof 
400?!^ ' 



t-- 



* See section {0) of Sir William Thompson's paper " On the Secular Cooling of the Earth." Thomp- 
son and Tail's Natural Philosophy, Vol. I. Part II. Appendix (D). The coefficient refers to the unit of 
bulk of the substance of the earth's crust. 
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or about 100,000,000,000 years. When applied to the case of the earth, there- 
fore, the series (10) and (11) are too slowly converging for the purposes of com- 
putation except when the time is very great. 

To overcome this practical difficulty and obtain an expression which will 
readily give the temperature u, of any point of the sphere at any time, we pro- 
ceed to transform the series (10.) 

4. In the transformation which follows, and again in the sequel, we shall 
make use of a certain Eulerian integral or Gamma function whose value may 
here be briefly demonstrated, the reader being referred for a fuller investigation 
to the better treatises on the integral calculus. The integral is 



e ax cos ax = - — e 

' 2(1 

o 

To prove the equality of the second and third members of this equation, differ- 
entiate with respect to /3 and then integrate by parts. Thus 

CO 

7y^ = — \ e xdx sm px 

o 

_ p-"'"'^sin/9^ 1 _ ^ r, - a^x^,i^ cos ^x. 
L 2a'^ J 2(^J 



That is 



or 



o o 

3v ^95/3 _ 

V 20?' 



whence log v =^ — -^ + log ^0. 

or J,^^^^-(i/3/«)^ 

To determine the constant v^ we observe that when /9 = o, 

00 

r —a'x'' J x/n 

V = v„= \ e ax ^ - — . 

" J 2(1 

o 

Therefore ^ 

v=- — e and e ^- = -^- 1^ axcosBx. 

2a y^TzJ 



OF INITIAL UNIFORM TEMPERATURE, ETC. 



79 



5. Applying the preceding integral to the «* term of equation (10), putting 
/32 = «V and 4a' = r,' j (a't), 



we have 



a-i/{7zt) J 
o 

Hence equation (10) may be written 



ru=. ~7T^ \e ^ ' dx Z ( — 1)"+'« 'cos«;rsm«7r — . (12) 

° n = \ 

Now, for brevity, put 
Then (12) becomes 



x^ylit and d^Ttrjr^. 



2rntln C —?-(, 2^2(402^2/)-! , V / \»+i -1 • a 

(13) 



_ ^oX r -r„V(4«^'rV)■ 



^ 



+ 2' (— i)" + '«-isin«0' + (?) 
« == I 

— 2" ( — !)" + '«"' sin «(_;' — d) 



Again, for brevity, put 



P = 



an'^/int) ' 






and denote the difference of the two trigonometric series in the right-hand mem- 
ber of (13) by Qx- Then, the expression to be evaluated is 



WW 

= P p-^'"" Qxdy. 



(14) 



6. We must now determine the value of Q^. For this purpose consider the 
well-known equation, 



W-g') 



= - -i- <?■ cos f -I- ^^ COS 2^ + g^ cos 3y -I- . . . 



I — 2g COS f -{■ g 2 

In this, for (f, substitute it ^ y ■\- f. Then we have 



W-g') 



2 \ ^ 2gco%{y -\- >f)-\r g 



-.^=gcos (j + f) — g"^ cos 2 (j/ -t- ^) + . . 
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Multiply this by d<p and integrate between the limits + d and — 6. The result is 

+ \g^ sin 3,(y + 0) — \g^ sin ^{y — d) 

- + 

= Q,sAy. (15) 

Comparing the series in (13) and (14) with that in (15), it appears that Qi^ Q 
when^= I. Hence equation (14) becomes 

-HO -0 O (,6^ 

The first term in the second member of this equation presents no difficulty. Its 
value is 

2p • 

But in the second term the element-function of the integral with respect to f is 
zero when g ^ i, except when cos [y -\- (f)-^ — i. In this case the element- 
function becomes 

I —g"" _ I +g _ 



{i—gf i—g' 



for g=i. 



The value of this integral with respect to <p is therefore the value of the single 
element — —^ dtp when ^ ^ i. To find its value the most direct process ap- 
pears to be the following* : — 

Replace cos {y -\- (p) by its equivalent i — 2 sin^ ii-f '^ f)- Then, after a 
little reduction, the integral becomes 



y, 1+ [i-^]\anH(j' + f) 



— rt - , . 

I +g 

Now every element of this integral is zero for ^ = I except that for which 

*The mathematical reader will recognize in the integral to be evaluated one which Poisson and other 
writers have used in proving the law of development of a function in a series of periodic terms. The dif- 
ficulty, if it be such, occurs in another form also in the proofs of that law given by Dirichlet and others. 
I am not certain that the evaluation given in the text will commend itself; but it seems to be more direct 
than, and quite as obvious as, Poisson's process. 
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J/ -{- (f = (2« + i)^-, n being any integer or zero; and the value of the excep- 
tional element is dependent on ti only so far as it determines within what limits 
jf can make J' + f = (2« + i) t:. We have, therefore, simply to consider what 
value the above expression has when y -{- f varies between {in + i) ;r — i and 
(2« + i)'^ + '' ^ being an infinitesimal. The indefinite integral is 



arc tan \j-Zf~. tan i ( J + <f)\ I 



I - — p- 
and since, when (i — g) is infinitesimal, — — - tan ■^■(j' + <f) varies from — oo 

to + 00 asj + f varies over its range, the proper limits of integration are — oo 
and + 00 . Hence we conclude that 

^ r (i-^Vy ^ = fare tan fL^tanK^ + ^l f = .. 

^J I +2^cos(/+ ^)+^' L li+^ 2 v^ 1 r/j J 

— —00 

It remains to find what values of / + f will give cos {j-\-<p)= — i, or 
what values of _j/ will make j + ^ an odd multiple of tt. Since (p may have any 
value between — and + 0, the range of values of jy, which makej + ^ an odd 
multiple of ff, will be from 

(2«+i);r — ^ to {2n-\- i)7r + d, 

n being any integer or zero. That \s,y must have values lying between 

n — and n -\- d, 
3n- — d and 3^ -f- d, 

t,n — d and 57r -\- 0, 



7. Recurring now to equation (16) and attending to the limits of j just de- 
rived we have 

/.„ = ^^^h^ — />;! f Ce-"'"' dy -f- Ce-"'^' dy + . . . 1 . 
Restoring the values of P, p, and d, and replacing py by z, there results 



i'\ — ^ 



= "'o-~"~ [ /^-" dz+ i^e-''dz+ . . . ] . (17) 

(18) 



2(7 1// 2ay t 

For brevity, let us put ;«„ = r„l{2a-\/t), 

in = rl{2a\/^i). 
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00 

Then observing that — 2— » i e~" dz =■ r^fi,^, 

o 

we get the following equivalent forms from (17) : — 

"'0 + >» S'^o + '« 



ru = ru, - ^ [f^"' ^^ + f^~" ^-^ + • • • 1 ' (19) 






m^—m 3m f, — I 



ru = ru,- r,u, [ 1 - A^je''^ dz — -^^e"' dz - . . . j , (20) 



\/nJ 1/7 

o nif, -\- m 

r K -u) = ?^ ^fe-' dz + fe-"dz+. . . ] , (21) 

(Wd — OT 3/«„ — m 

m„ — m yn^ — m 

r{u,-u) = r,u,\^i-^je-''dz-^je-'Wz- . . .J. (22) 

o /«(, 4- w 

Equations (19) and (20) give the product of the radius r and the excess of 
temperature u of any shell of infinitesimal thickness at any time. Equations (21) 
and (22) give the product of the radius and the fall in temperature of the same 
shell during the time /. 

8. It will now be well to determine whether equations (19) to (22) satisfy 
the conditions (i) to (5) and also to test their equivalence with (10). 

Differentiating equation (20) or (22) with respect to / we find 

"^ ,/;r ■ dt 



|/;r ' dt 



+ 



But, from (18), m^ — nt^-^ /-*, 



Hence 
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and ^^ =-^ [K - »«) ^-<"'>-""' - K + ;«)^ -<»» + »>' 

+ {im, -m)e- (»»» - ""' +...]. (23) 
Differentiating (20) with respect to r, there results, 

5r 






g (;«„ + ni) 

3r 



p-Cn^ + m)^ 



1 + 



By means of (18) 



3 (;«(, — m) I 

dr 2«]// • 



d{im^- 


m) 






I 


5r 








2ai/r 


5 (;«„ + 
3r 


m) 


= 


+ 


I 
2a^r 



Therefore 

9 (ru) 



^)?^L_ [-^-(».o -->'_(_ ^ -(»„ + »)' _|_^-(3».,- ».)=_|. _ _ .]-(24) 



5r "•*» «i/(;r/) 
Differentiating the last equation again, we find 



'■(ru) 



3r-^ 



a'ty TZ 






\- 



(25) 



Dividing (23) by (25), we get -p^ = «'-^. 

which is equation (i). The equations (19) to (24) therefore fulfill the first re 
quirement. 
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If we make r = r,,, equation (19) becomes 



r^u = r^u„ -^' I e 

o 



' fr« 



= ''o«o — 'V'o = O, 
and hence u = o. Condition (2) is thus satisfied. 
If we make r = 0, equation (20) gives 

ru = o — r^u^ y I — -i- J £-^' dz^ =0, 
o 
as required by condition (3). 

When t ^o, equation (20) becomes 



[^-v-J'-'^''] 



o 

or « = «o . 

which is in accordance with condition (4) . 

Finally, equations (19) to (22) should satisfy condition (5). But if we make 

^^ 00 in (21), the result is 

o o 

o o 

which is ambiguous. The same ambiguity appertains to equation (22). This 
difficulty will appear also when we seek to determine the temperature at the cen- 
tre of the sphere after an indefinitely great, or infinite, time. But it will be ob- 
served that if we can derive an expression from equations (19) to (22) for the 
temperature at the centre of the sphere, and show this expression to be zero 
when ^ = 00 , then those equations must satisfy (5). 
If we make r = o in (20), it becomes 

I — I 

u = u^ — rji.. 



o 
But the true value of this expression is readily found by differentiation to be 

ti = uJi —^{e-""' + e-"""'" + e-^'""' + . . • )]• (26) 

r = o. 

Now since m^^^^ rj 2a]/ 1, the last equation assumes the indeterminate form when 
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00 . But from the Eulerian integral of § 4 we have 

00 
- Ce-^^"'"''^'dx cos fix. 



-«v».o-' 



Making « successively equal to i, 2, 3, . . .in this expression, and substituting 
the equivalents in (26), we get 

00 

«< =: Kji [l — ^ I ^-(5^ '"■»'" ^;f (cos X + cos ^X + COS ^X -\- . . . )] . 
O 

By reference to the series used in § 6 it is seen that when ^ = i, 

tV s_J^ ^ s_j_ — cos ;i: + cos %x -f- cos 5;ir + . . . . 

I — 2^cos;r + ^'' I + 2^cos;ir + ^'' "^ ■' ' 

Hence we have when _^ = i , 

" L ?rJ i^i — 2^ cos JT + ^g-* I + 2^cos;r + ^''J J ^ '-^ 

o 

When / is infinite, m^ is infinitesimal, and we may for this value of / restrict the 

upper limit in (27) to an infinitesimal i, since for all greater values of ;f, ^-(J-"'""")' 

vanishes. In this case, also, the second term under the sign of integration in 

(27) may be dropped, by reason of the infinitesimal factor (i — g^). Therefore, 

when t^ CO and^= i, (27) becomes 

— i 

Since i may be assumed to be of a lower order than ;«„, the value of the factor 
^-(j.t/i'io)'' within the limits of integration is unity. The value of the integral 
then, as shown in § 6, is tz. Therefore we conclude that when r ^=. o and / ^ 00 , 

« = «o[i — (i /-)-] = o, 

and hence that equations (19) to (22) satisfy condition (5). 

9. Having established the complete equivalence of equations (19) to (22) 
and (10), we may now consider some practical features pertaining to their appli- 
cations. In the first place it will be observed that the definite integrals in either 
of equations (19) to (22) form an extremely converging series for all but very 
great values of the time. So converging are these series that, for the case of the 
earth, the first term will suffice when / is as great 100,000,000,000 years. When 
/ exceeds this time, recourse may be had to the original series (10), which then 
becomes highly converging also, as the following numerical examples will show. 

10. If we confine the time to relatively small values, it will suffice to use the 
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first term of (19) to (20). If, further, we restrict r to values which do not differ 
sensibly from r,,, we may write (22) in the form 



la-^/t 



«o — « = «« [ I — ~^S^~'^ ^^ ) • (^9) 



This agrees with an expression which Sir William Thomson has used for the 
fall of temperature, in his paper. On the Secular Cooling of the Earth* 

II. To illustrate the numerical application of the formulae derived, it will be 
well to assume such a value of the time that for the earth both sets of equations, 
i. e. (10) and (11), and (19) to (24), and (26) may be applied. 

As a first example let 

;- = \r^ and c?(Tt / r^t = i . 
The value of t required in the case of the earth to give this last equality is, as we 
have seen in §3, about one hundred thousand million years. With these data 
equation (10) gives 

u^^{e-^~\e-^+ . . .). 

log No. 

<?"' 9-5657055 — 10 +0.3678794 

e~^ 6.09135 — 10 I — O.0000411 

(e-'—^e-'-j-. . .) 9.5656570-10 0.3678383 

4 0.6020600 

^ 0.4971499 

«/?<„ 9.6705671 — 10 0.4683463; 

hence u = 0.4683463 Ug. 

*Thotnson and Tail's Natural Philosophy, Vol. I. Part II. Appendix (D). 

The law of cooling used by Thomson in this paper is that which applies to the diffusion of tempera- 
ture in an inHnite solid on the supposition that at the beginning of the time the temperature had two dif- 
ferent constant values on the two sides of an infinite plane. The formula which expresses this law is 



■2V c 

J' = i^o + ^ ,. I ' - " d^, 

o 
in which v^ is the half sum and Fthe half difference of the initial temperatures ; x is the distance of any 
point whose temperature is v, from the initial plane, the sign of x being positive for the hotter side and neg- 
ative for the cooler side of the plane ; / is the time from the initial epoch and '/. is the same as our a^. 
Now when /■==o the formula gives v = Vi^^ V, and hence the fall of temperature in any time / is 






This agrees with (29) if we make V= «„ and 
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Again, when ^'(iz J r^y( = i. 2ai/t^ zr^^jn; and when r= ^ra, equations 
(i8) give 

m = ^t:, 

m^ — in ^^s- ^ 0.7853982, 

m^-\- m =\-K= 2.3561945, 
3w„ — m ^^71 = 3.9269908, 
3^0 -\- m =\n= 5.4977871. 

s 

From the tables of the integral I e~'' dz given in Oppolzer's Lehrbuch zur Bahn- 

o 
besHnitnung, Vol. II. Table X, we find for the first two terms in (21) 

^e~"dz = 0.2355828.78, 

j e-"dz = 0.0000000.25. 

The sum of these integrals is 

log 
0.2355829 9.3721438 — 10, 
4. 0.6020600, 

i/tt 0.2485749, 

^5!^=^, 0.5316538 9.7256289—10. 

Hence by (21) «o — « = O.5316538 «o. 

or u ^ 0.4683462 Wo . 

This agrees with the value derived above from (10) within one unit in the seventh 
place of decimals. It will be observed that the contribution of the second term 
in (21) is less than one unit in the .seventh place. The first term of (21) would 
therefore suffice when the time is much greater than 100,000,000,000 years for 
all practical computations relative to the temperature of points within the earth 
12. As a second illustration, let us compute the temperature at the centre of 
the sphere when ?«o of (18) and (26) has the value 



zay't 
This gives t = -^ , 
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which for the case of the earth gives t = 273,000,000,000 years. We also have 
for use in (11), when Wo = i, 





UoJ ~4 




ence the computation by (i i) 


runs thu.s : — 






No. 


log 


1st term. 


\oge 


9.6377843 




_2 

4 


0.9942997 
0.6020600 


+ 0.0848050 


log.'i-)^ 


0.0300240 
1. 07 1 5785 
8.9284215 


2d term. 


log^(i^)' 
4 


0.03002 
0.60206 




log^-^^ 


0.63208 
4.2863 


— 0.00005 1 7 


e—' 


57137-.- 



10 



10. 



10. 

The third and higher terms are insignificant. The sum of the first two is 
0.0847533. Hence by (11) 

u = 0.1695066 «o- 
According to (26) the computation is as follows : — 

log 
+ 0.3678794 e-^ 9-5657055 — 10 

+ 0.0001234 e-"* 6.09135 —10 

+ 0.3680028 ^-' + ^-'+. . . 9.56585 1 1 — 10 

4 0.6020600 

V^ 0.2485749 

0.8304935 (4/|/7r)(^-' + ^-9_^ , . . ) 9.9193362 — 10. 
Therefore by (26) u ^ u^{\ — 0.8304935) 

= 0.1695065 Utj. 

Thus the results derived from (11) and (26) for the temperature at the centre of 
the sphere agree within the limits of precision attainable with 7-place logarithms. 



